124 


2.5A THE AR VS (AR> 2 QUESTION - THE PULSE-LENGTH DEPENDENCE 
OF SIGNAL POWER FOR FRESNEL SCATTER 

W. K. Hocking 

Max-Planck-Institut fur Aeronomie 
D-3411 Katlenburg-Lindau, FRG 


It has been proposed that the enhanced echoes from the atmosphere observed 
with a vertically pointing radar are due to reflections from horizontally 
stratified layers. The general case in which there are many closely spaced 
layers at random heights has been called "Fresnel scatter". The variation of 
received power with transmitter pulse length is examined for various models of 
Fresnel backscatter. It is shown that for the model most often used in previous 
work, the power is proportional to the pulse-length (Ar), and not to the pulse 
length squared. However, for more general models a pulse-length dependence more 
complex than either (Ar) or (A) 2 is found. 

1 . INTRODUCTION 

Radar backscatter at VHF from the troposphere and stratosphere shows at 
times evidence of weak partial reflections from extended horizontal 
irregularities. These irregularities are at least a Fresnel zone in horizontal 
extent, and fluctuate in the vertical by less than about A/8 over this 
horizontal distance* Here, A is the radar wavelength. This type of reflection 
is in addition to scatter due to turbulence-induced irregularities (e.g. , GAGE 
and GREEN, 1978; ROTTGER and LIU, 1978; ROTTGER, 1980a). GAGE et al. (1981a) 
have proposed that these scatterers occur at random heights in the atmosphere, 
and have then, using this simple assumption, proceeded to determine the expected 
dependence of backscattered power on the radar and atmospheric parameters. GAGE 
et al. (1981a) will be denoted by GBG here. The model was also discussed in 
GAGE and BALSLEY (1980), GREEN and GAGE (1980), GAGE et al. (1981b) and BALSLEY 
and GAGE (1981). The formula which GBG produced took the form 

2 2 

a M 9 9 

P„ ■= (Ar) (1) 

R 4X 2 r 

This formula is only relevant for the case in which the same array is used for 
both transmission and reception. P R is the received power, a is the array 
efficiency, is the peak transmitted power, A 0 is the array effective area, 

A is the radar wavelength, r is the range of the scatterers, M is the mean 
generalized refractive index gradient, and F(A) is a "calibration constant" 
which must be determined empirically for each radar. The term (Ar) represents 
the pulse width. 

Most of equation (1) is intuitively reasonable, but the (Ar) 2 terms 
appears to be odd. In this paper, the procedures adopted in obtaining this 
(Ar) 2 dependence will be carefully re-examined. It will be shown that there 
were errors in this original formulation, and that a proper treatment leads to 
a (Ar) dependence. 

This paper will primarily present the arguments for and against the (Ar) 2 
formula, although some mention will be made of generalizations of the Fresnel 
model. A more complete discussion has been presented elsewhere (HOCKING and 
ROTTGER, 1983). 
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2. PHYSICAL PICTURE 

A simplified view of the model presented in GBG is presented in Figure la* 
We will begin by discussing this simple model, and then will generalize it to 
gather complexity. 

Imagine that a square pulse of duration AT is transmitted upwards into the 
atmosphere, and at some time t 0 /2 the pulse is centred at a height (t Q is 
the time for the pulse centre to go to height z and be reflected back to the 
ground). Consider a vertical region of length cAT/2, centred on z 0 , and assume 
that within this volume, there are seven reflectors, of equal reflection 
coefficient but at random heights. Each reflector will reflect the pulse for a 
time duration AT, and at time t 0 = 2z 0 /c (where c is the speed of the radio 
waves), some part of the pulse will arrive back at the ground from each of these 
7 reflectors. No signal will arrive at time t 0 from reflectors outside of this 
region, (i.e., The received signal is a convolution between the pulse shape and 
the reflection coefficient profile.) The seven reflected signals will have 
approximately equal strengths, but because the reflectors have random heights, 
each signal will arrive back at the ground with random phase. The resultant 
signal may be described by the dark vector in the right-hand diagram of Figure 
la; that is, it is the sum of 7 vectors of equal strength but random phase. 

This is simply the classical two-dimensional random walk problem, as first 
described by RAYLEIGH (1894). It is well known that the modulus of the 
resultant vector of the two-dimensional random walk problem has a ,f Rayleigh 
distribution", and that the mean square length of the resultant vector is 
proportional to the number of contributing vectors. Thus if we double the 
pulse length to AT* and the mean number of reflectors per unit height remains 
the same at all heights, then we have approximately 14 reflectors in the new 
length CAT* / 2 in Figure la. As a result, an approximate doubling of the square 
of the resultant vector can be expected when the pulse length is doubled. 

Of course in the above discussion we dealt with small numbers of randomly 
phased vectors, and strictly speaking the Rayleigh distribution is only relevant 
for large numbers of vectors. Nevertheless, even for the cases of these small 
numbers of vectors, the mean power is still proportional to the number of 
vectors, provided that the reflectors are allowed to fluctuate vertically in 
time (so that each reflected component has a uniform phase distribution between 
0 and 2 tt c ) , and that the mean power is calculated over a long time interval. 
Naturally, however, the fluctuation in power about the true mean (relative to 
the true mean) will be smaller when larger numbers of reflectors contribute. 

More to the point, however, the above problem is only illustrative, and is 
unlikely to properly model the real atmosphere. Therefore, let us increase the 
complexity of the model. The above analysis at least gives one an intuitive 
feel that the power should be proportional to the pulse length. 

A more general model is represented by Figure lb. In this case, many 
reflectors are assumed to exist within one pulse length, but they are allowed to 
have varying reflection coefficients. This situation is analogous to that 
assumed in GBG. The situation is now far more complex than the classical 
random-walk problem. Nevertheless, HOCKING and ROTTGER (1983) showed that by 
dividing the reflectors into subsets of equal strength* it could be shown that 
the resultant vector will still be proportional to the number of contributing 
vectors, provided that the amplitude distribution of these vectors remains 
unchanged, and their phases are genuinely random. BECKMANN (1962) has 
considered the problem more rigorously and more generally, and has shown that 
the vector sum of a large set of vectors }, which have an arbritrary ampli- 
tude distribution but random phases distributed uniformly between 0 and 2^, is a 
vector with a Rayleigh distribution of amplitudes. Furthermore, BECKMANN (1962) 
has shown that the mean squared length of the resultant vector is proportional 
to the number of contributing vectors, in line with the above discussion. These 
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Figure 1. Pictorial description of reflection from a group of reflectors 
distributed randomly in height. The pulse is illustrated to the left in 
each figure, and the reflectors and their strengths are indicated by the 
horizontal lines. 

results are also consistent with NORTON et al. (1955). 

Therefore, it may be expected that the received power is proportional to 
the pulse length. 

The above arguments also apply if an arbitrary form of pulse shape is used, 
rather than a square pulse. Provided the pulse shape is sufficiently long that 
many reflectors contribute to any signal, the picture is still similar to Figure 
lb, but the amplitude of the pulse may change within the region cAT/2 . This 
simply weights the reflection coefficients, but the signals contributing to the 
total power at any instant are still due to signals reflected from a range of 
reflectors. These contributing signals are still uniformly distributed between 
0 and 2tt radians in phase, and the effect of the pulse is simply to modify the 
amplitude distribution of the component vectors. Thus the results outlined 
above still apply. 

We may now make a general statement. If we have a sequence of vectors 
{£. }, which have an arbitrary amplitude distribution and a uniform phase 
distribution (0-2ttc) , and this sequence is multipled by an envelope function E, 

then the vector sum of the resultant vectors {E. s . } obeys the relation 

— x i 

s 2 = W E , (2) 

where W E is the width of E, defined in any manner, and S 2 is the mean square 
vector sum. This relation is true for any specific envelope shape, but cannot 
of course be used to compare powers between different envelopes. 

It should be pointed out that if one or two of the specular reflectors are 
much stronger scatterers than all the others, the above statistical treatment is 
no longer valid. These cases require special consideration (e.g., RICE, 1944, 
1945; BECKMANN, 1962), but were not considered in the model of GBG and so will 
not be considered here. 

3. GBG TREATMENT 

In this section, the treatment adopted by GBG will be briefly outlined. 

For a more detailed treatment, the original paper could be consulted, as the 
description given here will be largely qualitative. Nevertheless, the principle 
is so simple that the "pictorial" treatment given here actually describes the 
model adequately. 

The approach adopted by GBG goes as follows. The reflection coefficient 
profile r(z) can be considered as the sum of many sinusoidal oscillations, of 
varying vertical scale, and varying amplitude. For example, the curves a, b and 
c in Figure 2 represent three of these. The amplitudes of these various scales 
can be found simply by Fourier transforming r(z). 




127 



Figure 2. Typical profiles of refractive index n(z), reflection 
coefficient r (z) , and 3 examples of Fourier components of r(z) . 

Then GBG state that the pulse is comprised of only one frequency, so only 
one of these vertical scales is important — namely, the scale with a node-to- 
node distance of A/2, A being the radar wavelength (i.e. , this is the Bragg 
backs cat ter scale). The received signal amplitude at the ground is proportional 
to the number of oscillations in length cAT/2. Doubling the pulse length 
effectively doubles the number of oscillations of this Bragg scale, increasing 
the amplitude received at the receiver by a factor of 2, and therefore the power 
by a factor of 4. A generalization of this discussion clearly suggests that the 
received power is proportional to the square of the pulse length. 

4. THE ERROR IN THE GBG ARGUMENT 

The argument in the previous section contains one error, and this is in the 
description of the transmitted pulse. GBG stated that a pulse consists of only 
one frequency, but by definition a single, pure frequency must be infinite in 
extent. A pulse comprises a spectrum of Fourier components, centred on the 
central frequency. As a result, a pulse comprises a range of wavelengths, and 
so in the description outlined in Figure 2, a finite spectrum of Fourier scales 
must produce backscatter. Since r(z) is a random function of height, the phases 
of these contributing Fourier components are random. Each scale therefore 
produces a reflected signal, and each signal arrives at the ground with 
different phase. These signals have random phase, so a "random-walk” type 
problem again results. 

It can be seen that a proper analysis is more complex than the simple 
description given by GBG. In the following section, the pulse- length dependence 
of the scattered power will be re-derived from the point of view of considera- 
tion of these various scales. It will be seen that the treatment given by GBG 
is inadequate, and the results of section 2 will be reinforced through this 
alternative approach. 

5. QUANTITATIVE TREATMENT 

Suppose that the pulse field strength at time t and height z is given by 

z~' 1 . g ( t-z / c ) . exp { jw ( t-z/ c )} , (3) 

where m is the carrier frequency, and g describes the pulse envelope. In this 
simple description, it has been assumed that the pulse travels at a speed 
c(= the speed of light in a vacuum), and absorption has been ignored. 
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For simplicity, the z" 1 dependence will be ignored. Let the pulse at z=0 
be written as 

2ja)g 

jwt c 

g t (t) e ■« g z (£)e * gp(^). (4) 


The function g t defines the pulse envelope, we have allowed gp to be in general 
complex, and £ = ct/2 is a length coordinate. We will consider only the case 
of g t symmetric about its maximum, as this is almost always valid for real 
experiments. The following results are probably true generally, independent on 
this symmetry refinement, but these asymmetric cases will be ignored for 
simplicity. Let us also associate a phase with r(z), where the phase is deter- 
mined by the height of the reflector above the ground, z. Then r(z) can be 
regarded as a complex profile, r(z). After backscatter from the reflection 
profile r(z), the signal received at time r 0 can be shown to be given approxi- 
mately by 


a(z o ) - r(z Q )* £ p (z Q ) 


“ / £(z) & (z o -z)dz, (5) 

—CO * 

where z q « cxjl. 

That is to say that the received signal is a convolution between x{z) and 
jjp ( z ) (e.g. , AUSTIN et al., 1969). It is convenient to work in the spatial 
domain, which is the reason that z Q has been used. The value z G can be approxi- 
mately regarded as the height from which most of the scattered signal received 
at time t g was reflected. 

Now introduce the functions A, R, and G, defined as the Fourier 
transforms of the functions a,, x and &p* That is, 

a. ( z ) ^ A ( 0 

r(z)«-»-R U) 

& (z)++G(0, ( 6 ) 

where ? is the reciprocal coordinate of z, (C plays the same role to z as 
frequency does to time; the C coordinate will be referred to as "reciprocal 
space*'.) Then, by the convolution theorem (e.g., BRACEWELL, 1978), 

ACO = G(0. RU). (7) 

Thus the signal strength received at the receiver can be found in the 
following way. First, find x(z)> and then find its Fourier transform R(C). 

Then find the Fourier transform of the pulse, G.(c). If R and G are 
multiplied, and then reverse Fourier-transformed, the signal amplitude aXz) 
can be found. This description is identical to the description given in section 
4, except that in this case we began by assuming a convolution in the spatial 
domain, whereas in section 4 we went directly to the reciprocal space domain. 
This shows that the treatment in section 2, and the discussion in sections 3 and 
4, are in fact different ways of viewing the same problem. We must now complete 
the analysis in the reciprocal space domain quantitatively, to show that it does 
in fact produce a pulse- length dependence for power. Given that GBG chose to 
work in the reciprocal-space domain, the following section gives the form of 
analysis which they should have adopted. 
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Since r.(z) is a random function of height, then R(z; ) is a random 
function of C . Therefore the function A(c ) is a random function with an 
envelope described by G(s). Schematic examples of K and G are shown in 
Figure 3, where gp(z) is taken as a Gaussian function, so G is a Gaussian 
function centred on C * 2/X . 

In any physical experiment, it is normal to , \aix ,f the central frequency 
down to 0 Hz, and for convenience we will do this in this theoretical considera- 
tion. This simply means that C “ 2/X is shifted to C * 0. Figure 4a shows an 
example of A( 0 after such a shift has been performed, and this function is 
denoted by A Q (C). Also shown (schematically only) in Figure 4 is the 
amplitude, |aj, and phase, ^ , of .a(£), which might typically result after 
A 0 (c) bas been reverse Fourier transformed. Notice that no large variations 
in | aj or <f> a can occur over distances of z of less than about one pulse length 
of gp(z). This is because G(C ) defines a limited frequency band of non-zero 
values A 0 (C ) , so no frequencies outside this band can occur in a( C). 

In any real situation, xjiz) will change as a function of time, and 
therefore so will a(z). The powers |a(z)j 2 at any height z may be averaged 
to produce a mean over some time interval T. This gives the mean power at 
height z. Since t(z) is random, there is no preferred" z value, and after 
sufficient averaging, | a | 2 will be a constant, independent of z. Therefore it 
is only necessary to look at one height, and for convenience we choose z = 0. 

By definition, 


a(z) - / A(c)e j2lrzC dC 

—CO 

so for z * 0 


( 8 ) 



Figure 3. Schematic illustration of R(C) , the Fourier transform 
of the reflection coefficient profile. (a) is the real part, 

(b) the imaginary component. Graph (c) shows the Fourier trans- 
form of the transmitted pulse. 
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0 ALTITUDE (Z) 


RECIPROCAL SPACE £ 

Figure 4, (a) The function A Q (c) (see text), (b) Typical amplitude 

and phase which might be recorded at any instant after reflection 
from the atmosphere, as a function of height, z. 


CO 

a(0) = / A (C)de. 


(9) 


In other words, a.( 0) can be regarded as the “complex area" under A^C). 

If we take d£ = Ac as a sufficiently small constant, 

N N 

U(0)| = | AC l A^CC ± )1 (= |AC l <G(C)R(C)|) (10) 

n=l n=l 

This amounts to simply vector ially summing a set of vectors {AcA oi }. It is 
clear that we are confronted with a very similar problem to that in section 2 — 
namely, we have a random sequence of vectors {ACR(Ci)}» which we multiply by 
some envelope G(Ci), and then we add to produce resultant* We wish to know 
how the modulus of the vector sum varies as we change the envelope width. The 
only difference compared to section 2 is that here the vectors are functions of 
reciprocal space, whilst in section 2 we were dealing with vectors which were 
functions of z. Clearly, then, the results in section 2 apply, and we see that 
if we hold the peak amplitude of G(c3 fixed, and define the ’Vidth' 1 of G(£ ) 
as W G , then the vector a(0) obeys the relation (2): i.e., 

|a.(0)|^ « 0-1) 

The width W may be defined in any way (e.g., half-power width, e" 1 width, etc), 
provided the definition is invariant for the chosen function. 


Equation (10) deals with the width of G( C) . It is now necessary to 
determine how changing the width of the Tx pulse g.p( c) affects G( 0 . Two 
results from Fourier transform theory are first necessarily. Firstly, the width 
of gp is inversely related to the width of G: i.e.. 



(12) 


and secondly. 
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S. (0) « / G(C)dS. (13) 

** —CO 

oo 

(conversely, G(0) ■« /& (z) dz). 

—CO P 

We are now in a position to examine the pulse- length dependence of the 
received power* As seen in equation (10), we have the following "random-walk" 
problem. We have a sequence of random vector { A£R( }, and we multiply them 
by an envelope {^G (Cj) }• We know that the transmitted pulse & p (z) may change 
in width but must maintain constant peak amplitude gp(G). Changing the width 
of gp(z) affects both the width and peak value of G(c). The width of G(c) is 
inversely proportional to the width of g (z) (by (12)), This fact, together 
with (13), means that the peak value of G(c) must be proportional to the width 
of gp(z) when gp(o) is held fixed. Thus both the width and peak value of G(c) 
change. If we consider the rescaling of the function G(c) and keep its width 
constant for now, we see that this simply increases all the vectors 
(A£R(Ci) G(Ci ) } by a factor proportional to W gp . This must therefore rescale 
the total power by Wgp 2 times. Now, we must only consider the effect of 
changing the width of the function G(c). Equation (2) can be applied here, so 
it is clear that changing the width of G changes the power proportionally to Wq. 


Combining the above effects, we have 


| a( 0) | « w 


gP 


(14) 


for the case of unchanging pulse peak power, and using (12), 

|a(z) | 2 «c w gp . (15) 

(We have already shown that |a(z) | 2 = |a(Q) l 2 for all z). 

This proves that the mean square received power is indeed proportional to the 
pulse width, even when viewed from the inverse space domain. 


6. COMPUTER SIMULATION AND GENERALIZATION OF ASSUMPTION 


Computer tests have been done to test equation (1), since that equation is 
crucial to all the arguments presented here. A Monte Carlo approach was adopted 
(e.g., SCHREIDER, 1967). The details of these tests will not be given here. It 
is suffice to say that equation (2) was completely verified by these numerical 
simulations. 


This Monte Carlo approach also allowed a generalization of the assumptions 
made by GBG. In the troposphere the mean reflectivity decreases approximately 
exponentially with height (e.g., BALSLEY and GAGE, 1981). Therefore the 
situation of a pulse incident on such a reflectivity structure has been 
investigated. In such circumstances, varying the pulse width will vary the form 
of the amplitude distribution of the reflected signals, and so the pulse- length 
dependence for scattered power is no longer simply proportional to (Ar), The 
details of this simulation can be found in HOCKING and ROTTGER, (1983), but the 
results are summarized here with Figure 5. Suppose that the RMS reflectivity as 
a function of height is <r 2 (z)> 1 / 2 , and that <r 2 (z)> 1 / 2 /z takes the form 
exp(-z/H). Suppose that a Gaussian pulse of half-power full width h is trans- 
mitted. Then the received backs cattered power is a function of h/H, and follows 
the form indicated in Figure 5. Clearly for h > 0.5 x H, the power is no 
longer simply proportional to the pulse width. 

7. DISCUSSION 


In any experiment to test the pulse-length dependence of backs cattered 
power, various precautions are necessary, or else misleading results can ensue. 
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h/H 

Figure 5. Plot of received power as a function of the ratio of the 
pulse width to the scale height of the reflection strengths, for 
the case of an exponential decay in reflection strength with height. 

Firstly, it is important that the receiver has sufficient frequency band- 
width to accommodate the pulse. For example, imagine the situation of trans- 
mitting a Gaussian pulse, and using a receiver matched to the pulse. (This is 
normally done, in order to optimize the signal-to-noise ratio.) If the trans- 
mitted pulse gp(z) is now made narrower, the peak value of G(c) falls 
proportionally. But if the receiver bandwidth is not widened to accommodate 
the wider range of frequencies, then equation (14) becomes 

|a.(0)| 2 cc [W(g p ) 2 ]. 1, (16) 

and the received power appears to be proportional to the square of the pulse 
width. I h fact, in any real investigations of this reflecting process, the 
effective pulse is not simply the transmitted pulse but rather that pulse 
convolved with the impulse response of the receiver. This last point is 
important, and care must be taken in petforming receiver matching. The receiver 
bandwidth must not be just equal to the bandwidth of the Fourier transform of 
the pulse, but considerably wider. For example, suppose that the transmitted 
pulse is described by g(t), and the Fourier transform of g(t) is G(.ca), w being 
the angular frequency. Let the receiver response be also G(u) — then the 
effective transmitted pulse is not g(t), but rather g(t)*g(t) — or a function 
roughly /2 times wider than the transmitted pulse. The receiver response 
should be flat over ail non-zero values of G(w) in order that the effective 
pulse is the same as the transmitted pulse. 

Secondly, if the receiver is matched to the transmitted pulse on all 
occasions, it is interesting to look at the signal-to-noise ratio. Doubling the 
pulse width doubles the received power — but if the noise is constant as a 
function of frequency over the bandwidth of the receiver, and the receiver band 
width is halved, to match the transmitter, then the received noise power also 
decreases by a factor of 2. For the case of VHF radars, the main noise is 
cosmic noise, and this can be regarded as constant over the band width of most 
VHF systems. Thus the signal-to-noise ratio is proportional to the square of 
the pulse length. It is important in any experimental test of the preceding 
theory to measure absolute power, and not signal-to-noise ratios. 
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HOCKING and ROTTGER (1983) presented a preliminary test of the above 
theory, using experimental data from the SOUSY radar. These results are 
su mm arized in Figures 6a-c, Figure 6a shows the experimentally observed power 
profile for a 150 m pulse after averaging over a period of 50 min. Figure 6b 
shows the profile which would have been observed had a pulse of length 1,5 km 
been used, with peak power equal to that of the 150 m pulse. A factor 
(1500/150) has been removed from this figure for ease of comparisons. This 
profile 6b was produced by computer manipulations; the details are discussed in 
HOCKING and ROTTGER (1983). 


(a) 


(b) 


6 MARCH 1981 



6 MARCH 1981 



Pm 


Figure 6. (a) Mean power as a function of height, recorded with the 

SOUSY radar on 6 March .1981, The noise has not been subtracted; 
the noise level was about 5-8 dB. (b) The resulting profile which 
would have resulted from using a pulse with a coarser resolution. 
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Figure 6c shows a comparison of the 2 profiles — and it can be seen that 
they are in approximate agreement, now that the factor (1500/150) (i.e, , the 
ratio of the pulse resolutions) has been removed from the low resolution pro- 
file, This is support for an approximately ( Ar) power dependence. However, a 
great many more experimental results are necessary to properly test the theory. 
Also, agreement is not perfect in Figure 6c; this is discussed further in 
HOCKING and ROTTGER (1983). 

Figure 6d shows the power as a function of height and time during this 
recording interval, after the mean power profile for the period has been 
removed. Notice the existence of certain stable, well-defined echoes. These 
are not consistent with the "Fresnel Scatter" model of GBG, and their existence 
must be borne in mind. This point is discussed further in HOCKING and ROTTGER 
(1983). The Fresnel scatter model may have relevance to the atmosphere, but it 
is not always applicable. 


6 MARCH 1981 



6 MARCH 1981 


(d) 



Figure 6. .(c) A composite of (a) and (b) . The solid line shows 

Figure 6(a) taken at steps of 1500 m. (d) Details of the echo 
strengths as a function of time during the period used to form 
the mean profile 6(a). Darker spots indicate greater intensity 
The mean profile over the period has been subtracted, so these 
plots are "residual signal strengths". 



135 


CONCLUSION 

The Fresnel scatter model by GAGE et al. (1981a) has been critically 
examined. It has been found that equation (1) is in error, and the (Ar) 2 part 
should simply read (Ar). Appropriate adjustment of F(A ) is also necessary. 

In the more general case of an exponential decay of <r(z ) 2 > l / 2 Jz with 
height, a more complex proportionality results, and this has been illustrated 
with a numerical Monte Carlo approach. 
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